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ABSTRACT 
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under  a convexity  assumption.  In  particular,  a theorem  of  Korovkin 
type  is  derived  for  a class  of  approximation  methods. 
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SIGNIFICANCE  AND  EXPLANATION 

About  ten  years  ago,  U.  Grenander  of  Brown  University  proposed  that  it  should 
be  possible  to  create  a theory  of  patterns  or  shapes.  Such  theory  would  have 
applicability  in  pattern  recognition  and  classification,  which  is,  for 
example,  of  great  importance  in  the  utilization  of  computers  in  biological 
and  medical  research. 

A particular  aspect  of  this  endeavor  is  the  approximation  of  shapes  in 
the  plane  or  in  three-space  by  simpler  shapes,  e.g.,  by  balls  or  ellipsoids 
or  other  simple  geometric  configurations  which  depend  only  on  a.  few  para- 
meters. As  an  outgrowth  of  this,  one  might  want  to  so  approximate  a shape, 
such  as  a tumor,  as  it  varies  in  time,  thus  producing  a simple  shape  which 
varies  with  time  or,  a set-valued  function  of  time. 

A set-valued  function  is  considered  here  to  be  a map  which  takes  [0,1] 

j 

into  the  compact  subsets  of  E . A definition  of  continuity  can  be  invoked, 
and  this  raises  the  question  of  approximation:  can  a given  set-valued 
function  be  approximated  uniformly  by  a "simpler"  one  (chosen  from  a given 
family) . Borrowing  well-understood  techniques  from  classical  approximation 
theory,  we  show  some  of  the  possibilities. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


APPROXIMATION  OF  CONVEX  SET-VALUED  FUNCTIONS 


Richard  A. 


Vitale 


1.  Introduction.  Let  K denote  the  collection  of  non-empty,  compact  subsets  of  IRd . 

With  the  introduction  of  the  Hausdorff  metric,  given  by 

(1.1)  h(K1,K2>  = inf (e  » o|ki  E k2  + eB'  K2  - Kl  + F'B^  ' 

K can  be  regarded  as  a complete,  separable,  and  locally  compact  metric  space.  Here  B 
is  the  closed  unit  ball  in  Kd  and  addition  and  scalar  multiplication  of  sets  are 
defined  in  the  usual  way. 

A set-valued  function  F is  a map  from  [0,1]  into  K.  Such  maps  (and  more 
general  versions)  arise  in  a variety  of  contexts,  including  optimal  control  theory, 
mathematical  economics,  and  probability  theory.  Analytical  investigations  have  followed 
several  lines,  including  the  construction  of  a differential-calculus  (see,  for  instance, 
Artstein  [2),  Aumann  [4],  and  Matheron  [8])  and  the  investigation  of  selections,  namely 
vector-valued  functions  f : [0,1]  -*■  IRd  such  f(t)  e F(t)  for  each  t (Wagner  [18] 
provides  an  extensive  survey  of  this  area) . 

Our  purpose  here  is  to  present  some  initial  investigations  into  the  possibilities 
of  an  approximation  theory  for  set-valued  functions.  We  take  our  lead  from  traditional 
notions  and  begin  by  posing  the  question,  is  it  possible  to  approximate  a given  F by 
a "simpler"  one?  More  concretely,  we  may  look  for  linear  approximants  of  the  form 

n 


(1.2) 


j=0 


* .K  . = * K 
3 3 0 0 


+ K 
n n 


where  the  K.  are  fixed  elements  of  K and  the  ^ . are  scalar  valued  maps  defined 
3 3 

on  [0,1].  A new  ingredient  in  this  classical  formulation  is  that  (1.2)  must  be  treated 
with  some  care  in  combining  terms.  Note  that  although  {0}  is  the  identity  for 
addition  of  sets,  i.e.. 


K + {0}  = K , 

generally  no  additive  inverse  exists  (one  can  easily  verify  that  K + K*  - {0}  cannot 
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be  solved  for  K1  unless  K reduces  to  a point) . Moreover,  the  distributive  law 

iK  + 8K  = (ct+B)  K 

generally  fails  to  hold  (consider,  for  instance,  the  case  when  K » (0,ll  JR^) . It  is 
true  that  a restricted  version  of  11.2)  holds  for  convex  K,  namely 

(1.3)  aK  + 6K  = (a+8) K for  a, 8 > 0 . 

This  suggests  that  the  class  of  convex-valued  F may  be  an  appropriate  place  in  which 
to  begin  considering  approximation,  and  we  will  devote  our  discussion  to  this  case. 

An  outline  of  the  development  is  as  follows.  In  section  2,  we  present  notation  and 
generally  well-known  preliminaries.  We  take  up  Bernstein  approximation  in  section  3 to 
show  the  possibility  of  uniform  approximation  by  linear  approximants  of  polynomial  type. 
We  then  make  a brief  excursion  into  the  non-convex  case.  Section  4 presents  our  main 

result  - of  Korovkin  type  - on  positive,  linear  operators.  In  section  5,  we  return  to 

Bernstein  approximation  to  examine  some  of  its  other  features. 


I 


2.  K^.  We  denote  by  the  collection  of  elements  of  IK  which  are  also  convex. 

We  summarize  some  properties  of  K which  can  be  found  in  standard  references  (see,  for 

c 

instance,  Eqqleston  [7],  Rockafellar  (10],  and  Valentine  [14]). 

is  closed  under  addition  and  scalar  multiplication  of  sets  and  enjoys  the 
distributive  property  (1.3).  inherits  its  metric  from  K as  a closed,  separable, 

and  locally  compact  subspace.  Given  an  element  K,  we  may  form  its  convex  hull,  which 
is  in  K^.  The  map  K ■*  con  K is  continuous  and  satisfies  additionally 

contaK^+SK^)  = a con  + 8 con 

for  a , 8 >_  0. 

To  each  K e K__  is  associated  its  support  function,  given  by 

(2.1)  s(p,K)  = max{p-k|k  e K)  p e lRd,  |jp||  =■  1 . 

One  may  consider  the  support  function  to  give  a convenient  parameterization  of  the  family 
of  supporting  hyper planes  to  a set.  A set  K e and  a point  not  in  K can  always 

be  separated  by  some  hyperplane,  and  this  leads  to  the  useful  equivalence 

(2.2)  Kx  C K2  <=>  s (p,K1)  <s(p,K2)  Vp 

and  consequent  uniqueness  of  support  functions 

Kx  = K2  <=>  slp,^)  = s(p,K2)  Vp  . 

As  a function  of  p,  s(p,K)  is  continuous;  indeed  the  Schwarz  inequality,  together 
with  (2.1),  yields  the  uniform  bound  |s(p2>K)  - sfp^.K)!  <_  || p2  - p^  ||  ||k||.  Here  we  have 
used  the  symbol  [[ K [|  to  denote  the  norm  of  K which  is  equal  to  max{||k|||k  e K)  and, 
equivalently,  d((o),K). 

Evidently  we  may  use  the  map  K ►*  s(-,K)  to  embed  K in  the  Banach  space  B, 

c a 

of  continuous  functions  defined  on  the  surface  of  the  unit  ]Rd  ball.  Important  structure 
is  preserved  under  this  mapping: 


(2. 

.3) 

s ( - ,oK) 

= (IS  ( ' 

■ ,K)  a > 0 

(2. 

.4) 

s(-,k1+k2) 

- s ( • , 

,Kj)  + s(-,K2) 

(2. 

.5) 

h(K1,K2) 

■ llsl 

- s ||  (uniform  norm) 

<11 'll 

* II  S ( • 

■ ,K)  || ) . 

-3- 


Let  us  indicate  briefly  how  (2.5)  comes  about:  The  support  function  of  B is 
1 so  that  (2.3)  and  (2.4)  imply  s(p,l<2+FB)  = s(p,K2)  + e.  Together  with  (2.2) 
yields  k C k2  + eb  iff 

stp.Kj)  <_  s(p,K2>  + E for  all  p . 

The  analogous  expression  holds  iff  — K1  + For  inclusions  to  hold, 

(2.6)  Isfp.Kj)  - s(p,K2)|  < e for  all  p . 

The  infimum  of  all  e > 0 satisfying  (2.6)  is  at  once  h(Kj,K2>  and  ||  s^  - s 

(1.1)).  Taking  in  particular  K2  = {0}  yields  || K ||  = ||s(-,K)||. 

<T  [ IK]  and  <E[K  ] will  denote  the  spaces  of  continuous  maps  into  K and 
c 

tively.  Given  a map  F e I (IK)  we  denote  the  norm  of  F by 

H (F)  = sup{  1 1 F (t)  ||  1 1 e [0,1]  } 
and  define  the  related  metric  by 

H(F,G)  = sup{h(F(t),G(t))  |t  e [0,1]}  . 


identically 

this 


we  must  have 

, (see 

K respec- 
c 
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3.  Bernstein  Approximation.  Given  a set-valued  function  F,  we  define  its  n 
Bernstein  approximant  to  be 


th 


B (F;t>  = l b.  ( t ) F ( j/n)  , b.  (t)  = ("It^  (l-t)n"j  . 

n 3n  3n  I3/ 


It  is  straightforward  to  show  that  this  map  necessarily  lies  in  C[1K)  and,  indeed,  in 


d[K  ] if  Ft  (T  [ K ) . 

c c 


Theorem  1 . Let  F c ff[K  ).  Then,  as  n -»  “,  B ( F ; ) converges  uniformly  to  F (i.e. 


H (F  , B (F;  ) ) ->  0)  . 
n 


Proof.  We  use  the  Banach  space  embedding.  F e I[K  ] is  equivalent  to  the  continuity  of 


the  map  from  [0,11  into  B given  by  t s(*,F(t)).  Likewise,  a Bernstein  approximant 


of  F corresponds  to  the  map  t »■  l b . (t) s ( • ,F ( j/n) . Hence,  it  is  enough  to  show 


3=0 


the  uniform  convergence  (in  of  the  latter  maps  to  t h-  s(-,F(t)).  This  follows 


directly  from  classical  arguments  (see,  for  example,  Davis  [5]  transposed  to  a Banach  space 
setting) . O 

In  section  4,  we  shall  view  this  result  from  a more  general  perspective.  For  now,  let 
us  turn  to  the  case  when  F e d[K]  does  not  necessarily  have  convex  values.  Of  course. 


this  does  not  preclude  forming  B^(F;)  and,  indeed,  as  we  shall  see,  Bernstein  approxima- 


tion asymptotically  "convexif ies"  F. 


Let  us  digress  for  a moment  to  consider  a simple  example.  If  K = {0,1}  C it,  then 


1 n 

~ I K = i [K  + K + + K]  = {0,-,-,  . . .,1} 


3=0 


and  hence  hi  'j  K,  con  K — [0,11  ) 0 as  n -*■  «°.  This  "filling  in"  of  values  is 

j=0  ' 

typical  of  what  happens  when  non-convex  sets  are  summed.  The  following  result  quantifies 


this  behavior . 


Proposition  (Shapley-Folkman;  see  Arrow  and  Hahn  [1,  p.  396)).  Let  f K , j = 0,1, ...,n. 


be  such  that  ||K^[|£_M.  Then 


(3.1) 


l K . , l con  K . < Md 


3=0 


j=0 
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We  use  this  result  to  investigate  the  non-convex  case. 

Theorem  2.  Let  F e <E[K].  Then  in  any  subinterval  (f,l-fj,  0 < e < , B (F ; ) 

2 n 

converges  uniformly  to  con  F. 

Proof.  With 

n 

B (F;t)  = l b.  ( t)  F ( j/n)  , 
n . ^ in 

3=0 

we  identify  K.  = b.  (t)F(j/n)  in  (3.1).  Now 
3 jn 

|| K,.  II  < II F ( j/n)  II  |bjn(t)  I < H(F)  supfb^lt)  |e  < t < 1 - e.  j = 0,1 n} 


The  indicated  supremum  can  be  shown  to  be 


o(n"2). 


so  that  by  the  proposition 


h (B  (F;  t)  , B (con  F 
n n 


I -i|  1 

;t)  ) < H (F)  01  n 2 )d2  . 


Theorem  1 applied  to  B (con  F;)  and  the  triangle  inequality  yield  the  assertion.  □ 
n 

We  remark  that  the  result  cannot  be  extended  to  the  full  interval  since  at  each  end- 

( ‘I I 

point,  t = 0,1,  B (F;t)  = F(t)  independent  of  n.  Moreover,  the  01 n / bound  breaks 
n 

down  at  the  end-points. 


The  convexif ication  of  F by  Bernstein  approximation  is  undoubtedly  related  to 
theories  of  integration  of  set-valued  functions,  which  invariably  yield  integrals  with 
convex  values.  It  would  be  of  interest  to  make  this  statement  more  precise  via  a general 
investigation  of  the  behavior  of  linear  operators  on  set-valued  functions.  We  shall 
not  consider  this  problem  here  but  instead  present  another  example  which  shows  the 
difficulty  of  formulating  approximation  methods  in  the  non-convex  case. 


Let  F (t)  5 {0,1}  be  approximated  by  a piece-wise  linear  scheme 


F (t)  = ( [nt J+l-nt)  F + (nt  - IntJ ) F 11  nt-^  0 < t < 1 

n In  n — 


- F (1)  , t = 1 


Here  F (0)  = F (1)  = {0,1 },  whereas  the  sequence 
n n 

F (t)  = {0,nt  - Inti,  Int)  + 1 - nt,l) 
n 


even  fails  to  converge  for  any  t £ (0,1). 
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4.  A Korovkin  Result.  Using  (2.2)  and  the  positivity  of  Bernstein  approximation  of 

real -valued  functions,  we  see  that  for  F,G  e (T[K  3 

c 

F C G (that  is,  F(t)  C G(t)  , Vt)  =>  B (F;  ) C B (G;)  Vn  . 

— — n — n 

As  in  the  real  case,  this  suggests  that  a wider  class  of  approximation  methods  may  possess 

similar  convergence  properties.  Let  us  agree  to  call  a map  T : (C  ( IK  ] -+•  <E  ( K 1 

c c 

K -linear  if  T{aF+0G)  = aTF  + BTG,  a,  B > 0 and  K -positive  if  F C G *>  TF  C TG 
c — c — - 

(note  that  the  definition  of  linearity  is  restricted  to  the  natural  form  for  convex  sets) . 
We  then  have  the  following  Korovkin  result  for  such  maps. 

Theorem  3.  Let  {T^}  be  a sequence  of  -linear,  IK  -positive  maps.  In  order  that 
-*■  F for  each  F e it  is  necessary  and  sufficient  that 

(i)  T FCl)  1-  F(l)  i = 0,1,2  where  F(l>  (t)  = tXB 

n 

and 


(ii)  sup{H(T^F,F)  |F(t)  = K,  ||  K ||  = 1 } 0 . 

Let  us  remark  briefly  on  the  hypotheses  of  the  theorem  before  proceeding.  Condition  (i) 

is  reminiscent  of  the  vector-valued  formulation  and  is  perhaps  even  more  striking  here 

in  that  only  a fixed  shape  (i.e.  B,  the  closed  unit  ball)  is  involved.  Condition  (ii) 

asserts  that  the  T behave  uniformly  well  when  applied  to  "constants"  (including  the 
n 

case  F = F^  from  (i)). 

Necessity  of  the  conditions  is  the  more  direct  implication  to  verify:  (i)  is 

obvious.  As  for  (ii) , suppose  the  contrary.  Then  there  is  an  r > 0 and  a sequence  of 

K such  that  H(T  K ,K  ) > f (here  we  have  abused  notation  slightly  to  let  K 
n . n . n . n . — n . 

3 3 3 3 3 

stand  for  F where  F (t)  = K ) . Local  compactness  of  K and  the  uniform 
n . n . n . c 

3 3 3 

normalization  ||K  ||  = 1 assure  the  existence  of  a convergent  subsequence  of  the  K 

"j  nj 

Without  loss  of  generality,  suppose  that  K K . Then  by  the  triangle  ineauality, 

n «> 

H (T  K , K ) < H (T  K ,T  K ) + H (T  K , K ) + H (K  ,K  ) . 

nnn  — nnn«  n °o  oo  <»n 


Now  r = H(K  ,K  ) -*■  0.  Moreover,  the  twin  inclusions 
n oo  n 

K c K + e B 
n — 00  n 

K C K + r B 
n — n n 
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together  with  the  properties  of  imply 


TK  C T K + fTB 
n n — n “ n n 


TK  CTK  + e T B 
n — n n n n 

so  that  H(T  K ,T  K ) < e H(T  B)  -*■  0.  Hence  lim  H(T  K ,K  ) > £,  but  this 
n n n — n n n <■>  °°  — 

violates  our  assumption. 

The  proof  of  sufficiency  is  more  involved  and  will  require  some  preparation.  We 
begin  by  formulating  a quantitative  Korovkin  result  for  families  of  real-valued  func- 
tions. This  will  then  be  adapted  to  our  needs  by  invoking  the  Banach  space  embedding. 

Let  P be  an  indexing  set  and  let  E denote  the  collection  of  all  o = (o^),  P £ P- 
Here  we  have  denoted  by  < an  equicontinuous  family  of  real-valued  functions 

defined  on  [0,1].  That  is,  given  o = <0  >,  each  a £ C[0,1]  and 

P P 

(i)  3M  such  that  sup  II o II  < M < °°  , 
o p"  — a 

P 

(ii)  The  modulus  of  continuity  oo  (6)  = sup  sup|o  (t)  -a  (x)  | 

° It-* I <6  P P P 


satisfies  oi  (0+)  = 0. 
o 

E is  a normed  linear  space  under  the  definitions 


aa(1)  ♦ Ba(2)  = <ao(1)  + Ba(2>> 


|1  a ||  = sup  sup|o  (t)  | . 
t P P 

Moreover,  we  can  define  a partial  ordering  by 

a(1)  < a(2)  <=>  a(1)  (t)  < a(2)  (t)  Vp  c p,  Vt  e [0,1]  . 

P - P 

Now  let  us  consider  a subspace  Eg  C E and  a map  L : Eg  -*■  E.  We  say  that  L is 

..  (2) . (1)  (2)  . . . (1)  J (2)  r (1)  . (2) 

linear  if  L(cto  +Ba  ) = aLo  + BLa  and  positive  if  a -s  a =>  Ln  Lo 

For  convenience,  we  call  Eg  full  if  the  following  conditions  hold 


(i)  For  i = 0,1,2,  .o  e EQ  where  ^o^ft)  = t Vp  f P,  Vt  f (0,1]  (note  that, 

since  E is  a subspace,  this  implies  that  , .o  e E where 
0 lx]  0 

, ,o  (t)  = (t-x)2  Vp  £ P,  Vt  £ [0,1],  xe  [0,1]  fixed)  . 

I X j P 
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% 


(ii)  If  a = <o  > e E„,  then  for  each  fixed  x e (0,1), 

P 0 

,o  e I where  , .0  it)  = o (x)  Vp  e P,  Vt  e [0,11 
(x)  0 (x)  p p 


(loosely,  must  contain  enough  constants) . 


Finally,  we  define 
(4.1) 


y(o,L)  = sup  sup  [L  o]  (t)  - o (t) 
t p I (t)  P P 


We  are  now  prepared  to  state  a uniform  bound. 

C 

0 


Proposition.  Let  be  full  and  let  L : En  -*•  E be  a positive,  linear  map.  Then  for 


each  a e E 


||a  - La ||  1 (u)  [ II L0a  II  + 1]  + ylo,L) 

2 

where  u = sup  sup|[L  ,ol  (x) | . 

[x]  p 

x P 

Proof.  We  follow  an  argument  of  Shisha  and  Mond  [13],  who  have  developed  a 
similar  quantitative  estimate  in  the  case  where  P has  a single  element. 
Fix  o = (a  > e E . Then,  for  each  pep  and  6 > 0, 

p 0 


o (t)  - a (x)  | < id  (6) 
P P — a 


1 + 


(t-x) 


21 


Consider  one  of  the  two  associated  inequalities,  for  example, 

, 2' 


a (t)  < on  (6) 
p - a 


1 + 


(t-x) ' 


+ a (x) 
P 


Regarding  x as  fixed,  we  see  that  this  is  equivalent  to 

1 


a (6) 

^ o 


a + — a 
0 52  [X] 


+ a 
(x) 


and  hence 


Lo  ^ (1)^  ( 6 ) 


L 0°  + 52  L [x]° 
th 


+ L , 0 . 

(x) 


The  opposite  ordering  is  similar.  We  take  p components  evaluated  at  x and 
combine  the  two  resulting  inequalities  to  get 


By  assumption, 


[Lol  (x)  - [L  a]  (x)  | < u>  (6) 
p (x)  p 1 — o 
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(L0olp(x)  + (xl«lp<*>(  - 

6 


[L  o]  (x)  - a (x)  < Y(o,L) 

(x)  p p — 

and  two  applications  of  the  tria  igle  inequality  yield 


(Lai  (x)  - a (x) 
P P 


00)p<x)l  + JJ  fxl°]p(x)l 
6 


* y(o/L) 


Taking  sup  sup  on  each  side  yields 


- <V6>  II L 0CT 1 1 + \ n2 

6 


+ Y (o,L)  . 


If  P > 0,  we  take  6 = p and  are  done.  If  p = 0,  then  a limiting  argument  (see  [13]) 
similarly  yields  the  assertion.  □ 

Corollary.  Let  Z^  be  full  and,  for  each  n = 1,2, let  L^  be  a positive,  linear 

map  taking  Z into  Z.  If  L .a  -*■  .a  for  i = 0,1,2,  then,  for  each  a e r , 
u nil  0 

Y(o,L  ) •*  0 implies  L a ■*  a. 
n n 

Proof . In  view  of  the  proposition,  we  only  have  to  show  that  p2  = sup  sup|[L  a]  (x)j-*-0. 

n x p n x p 

2 2 

Note  that  each  component  of  [x]°  is  t - 2xt  + x (here  t is  the  free  variable) 
or  equivalently 

2 

a = a — 2x  a + x a . 

[x]  2 ZX1  X 0 

We  apply  L , take  pth  components,  and  evaluate  at  x: 
n 

[L  o]  (x)  = [L  a)  (x)  - 2x[L  ,a]  (x)  + x2[L  a]  (x)  . 
n ( x ] p n 2 p nip  nOp 

2 

Adding  and  subtracting  2x  appropriately  on  the  right  and  taking  absolute  values  yields 
the  bound 

(4.2)  |[Ln  (x]a)p(x)|  I |[Ln  20]p(x)  - X2|  + 2|[Ln  ^Ipfx)  - x| 

+ l[Ln  00lp(x)  - • 

Operating  with  sup  sup  on  each  side  of  (4.2)  then  gives 
x p 

% - II  Ln  2°  " 2°  I'  + 2 II  Ln  1°  “ 1°H  + HLn  0°  " oaH 
and,  by  assumption,  each  of  the  three  terms  on  the  right  tends  to  zero.  D 


We  now  adapt  these  considerations  to  the  proof  of  sufficiency  in  theorem  3.  Recall- 
ing the  identification  between  a set  K e ik  and  its  support  function  s(p,K), 

c 
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aF  + 8G  *-*  as(F)  + 6s  (G)  a,  6 >.  0 . 

A Kc -linear  operator  T : -*•  IfK^l  induces  a natural  map  L : C -+•  C via 


Ls  (F)  5 s (TF) 


which  obviously  satisfies 


L [as  (F)  + 6s  (G)  1 = aLs(F)  + 6Ls(G)  a, 6 >_  0 . 

In  order  to  apply  the  proposition  and  corollary,  we  need  to  extend  the  domain  of  L to 
a subspace  of  £.  Accordingly,  let  £Q  be  the  span  of  C (i.e.  all  finite  linear 
combinations  of  the  form  I a^s(F^))  and  define,  for  any  s(F), 

L[-s (F) ] = -Ls (F)  (=  -s(TF))  . 

With  this  done,  it  is  straightforward  to  verify  that  L : £ -*  E is  linear.  Moreover, 

since  F C G <=>  o(F)  ^o(G),  it  follows  directly  that  if  T is  K -positive,  then  L 
is  positive. 


Next  we  show  that  £^  is  full.  Obviously,  F 


s<F  > = i°  f V 1 = °'1'2- 


Further,  , .a  is  of  the  form  , .a  (t)  = E a . s (p,F . (x) ) . But  each  s(*,F.(x)), 

(x)  (x)  p 11  1 

regarded  as  a function  of  t,  corresponds  to  a constant  set-valued  function  F(t)  F^  (x)  . 

Hence  each  , .at  £„. 

(x)  0 


we  see  that  an  F e i[k  ] can  similarly  be  identified  with  its  family  of  support 


functions 


s (F)  = < s(p,F(.)  > 


where  the  indexing  set  is  P = {p|||p||=l}.  Now 

sup|  s (p,F  (t) ) | = 1 1 F ( t ) ||  <^H(F)  < oo 

P 


sup  sup) s (p,F (t) ) - s(p,F(x))|  = u (6) 

! t-x | <6  p 

where  w is  the  modulus  of  continuity  of  F defined  in  the  obvious  way.  Hence  s(F)  6 T. . 
In  fact,  the  collection  of  all  s(F)  form  a positive,  convex  cone  C in  E by  virtue 
of  the  identification 
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IK^  -linearity  and  -positivity  obviously  hold.  Moreover,  for  i = 0,1,2, 

B^tF^st)  = B^(tSt)-B  which  establishes  convergence  for  the  Finally,  given 

any  constant  F(t)  2 K,  Bn<F;t)  = F(t)  so  that  the  derived  y in  each  case  is  zero. 

It  is  equally  straightforward  to  establish  convergence  of  the  piece-wise  linear 
scheme  (3.1) . 


Bernstein  approximation  which  complement  the  uniform  convergence  result.  Some  have 


been  alluded  to  before  and  are  true  for  similar  approximation  schemes. 

We  begin  with  some  properties  which  follow  directly  from  the  support  function 
embedding  and  properties  of  Bernstein  approximation  in  the  real-valued  case. 

Proposition. 

(i)  K,  c F (t)  c K , Vt  =>  K c B (F;t)  C K,  Vt . 

1—  — 2 1 — n — 2 

In  particular,  n F(t)  CBn(F;t)  C conj’U  F(t)J  Vt. 

(ii)  F (s)  C P)F(t)  Vs  < t =>  B (F;s)  C p)B  (F;t)  Vs  < t . 

— — — n — — n — 

(iii)  f|S  * tj  C (D)  | [F  (s)  + F (t)  ] Vs , t 

=>  B [f;  — | c P)  i [B  (F;  S)  + B (F;  t)  ] Vs , t . 

n|  2 I — — 2 n n 

Property  (i)  is,  of  course,  a special  instance  of  the  positivity  property.  As  we 
have  seen,  this  is  a natural  extension  of  the  real-valued  case.  One  might,  however, 
try  to  argue  another  type  of  extension.  If  f (t)  > g(t)  Vt  then  the  Bernstein 
approximants  of  these  functions  share  the  same  ordering.  Alternatively,  one  could  say 
that  non-intersection  of  graphs  is  preserved.  Accordingly,  in  the  set-valued  case,  we 
might  ask  whether  non-intersection  - F(t)  l"l  G(t)  = 0 Vt  - is  maintained  for  approximants. 
The  following  example,  however,  shows  that  this  is  not  generally  the  case.  In  the 
complex  plane,  let  F(t)  = {e2l,lt}  and  let  G = {z|||z||  <_  e}.  Then,  for  each  t,  B^fFjt) 
is  a point,  which  for  t^  0,1  is  of  modulus  less  than  1.  Hence  if  e is  sufficiently 
close  to  (but  smaller  than)  unity,  F(t)  n G(t)  = 0 Vt  whereas  this  property  fails  for 
the  approximants.  It  is  possible,  however,  to  show  that  non-intersection  is  ultimately 
preserved  in  general. 

Proposition.  Let  F(t)  n G(t)  = 0 Vt.  Then,  for  n sufficiently  large, 

B (F ; t)  n B (G;t)  = 0 Vt. 
n n 

Proof . Let  r = inf  inf { ||  f -q||  I f e F(t),g  e G(t)}.  Compactness  and  continuity  ensure 
t 

that  e is  strictly  positive.  The  assertion  then  holds  for  n such  that 

e/2  > max { H ( F ,B  (F  ;)),  H (G, B (G ; ) ) } . □ 

n n n n 


J 
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We  turn  now  to  the  behavior  of  approximants  when  juxtaposed  with  mappings  of  the 

"background  space"  . If  M is  a d >'  d matrix,  we  can  define  a map  taking 

into  K by  K " UK  = (Mklk  6 kl.  The  following  easy  result  is  typical, 
c 

Proposition.  B (MF;t)  = MB  (F;t). 

In  particular,  Bernstein  approximation  commutes  with  projections.  Alternatively,  let  us 

consider  a continuous  one-parameter  family  of  matrices  M^,  0 < t < 1 (continuity  can  be 

assumed  in  any  reasonable  sense,  e.g.  in  the  Euclidean  norm).  Then  F(t)  = MtK,  for  a 

(one  might  think,  for  instance,  of  a continuous 

rotation  of  a fixed  figure).  As  well  as  uniform  convergence,  we  have  the  following. 

Proposition.  B (M  ;t)K  C B (M  K;t)  . 

c n t - n t 

Proof.  As  a consequence  of  the  general  inclusion  (M^+M^)k  CM^K  + M^K  (matrix  addition 
on  the  left,  set  addition  on  the  right) , we  have 


fixed  K e X is  an  element  of  I[K 
c c 


B (M  ; t) K 
n t 


y b.  (t)M . , KC  Y b.  (t)M . , K = B (MK;t) 
j=0  3”  3/«/  ~ m 3/n  n t 


□ 


Note  incidentally  that  B (M  ; )K  also  converges  uniformly  to  M.K  which  suggests 

n t 

further  comparisons  with  the  convergence  of  the  Bernstein  approximants. 

An  area  of  particular  interest  is  the  behavior  of  geometric  functionals  under 

Bernstein  approximation.  For  instance,  given  any  functional  : Kc  ■+  1R1  which  satisfies 

(ctK,  +8K_ ) = cup  (K, ) + &P(K„),  C1,B  > 0,  we  have  the  obvious  relation  <P(B  (F ; t ) ) = B (<p°F;t) 
1 2 1 2 — n n 

Examples  are,  for  fixed  p,  ^ (K)  = s(p,K),  the  extent  of  K in  the  direction  p and 

<f  (K)  = s(p,K)  + s(-p,K),  the  width  of  K in  the  direction  p.  In  the  plane, 

>p(K)  = per(K)  = perimeter  of  K is  another  example.  Here  a convenient  parameterization 

takes  p = (cos  0,  sin  0)  so  that  the  support  function  may  be  regarded  as  a function  of 

2tt 

the  angle  0.  Then  per(K)  = / s(0,K)d0  (see,  e.g.  [14])  and 

0 

n 2m 

per  (B  (F;  t) ) =5"  b.  (t)  / s (0 ,F  ( j/n) ) d0  . 

j=0  3"  JQ 

Nonlinear  functionals  naturally  require  individual  attention.  Occasionally  classical 


considerations  can  be  invoked,  as  in  the  following  bound  for  the  volume  of  B^fF.-t), 


-14- 


which  is  a straiqhtforward  consequence  of  the  Brunn-Minkowski  inequality  (see,  for 
instance,  [7] ) . 


Proposition,  vol  B (F;t)  » (B  ( (vol  FJ^'Stll^. 

n — n 

In  the  plane  an  explicit  expression  can  be  displayed  for  the  area  functional, 
sufficient  smoothness  of  the  support  functions 

n n 2n 

area  B (F;t)  = - £ £ b.  (t)b  (t)  / s (0  ,F  ( j/n) ) -r  (6  ,F  (k/n)  )d6 

" j=0  k=0  ]n  kn  0 


where  r(0,K) 


— — + I s(0,K)  (see  the  discussion  in  (11)  on  mixed  areas). 

36 


With 
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6.  Notes. 


§1.  The  present  study  was  motivated  in  part  by  earlier  work  of  the  author  and 
colleagues  in  related  areas  - in  particular,  approximation  of  plane  convex  sets  (16), 

[9]),  random  sets  ([31,  [171),  computational  considerations  ([15)),  and  modelling  of 
tumor  growth  ([16)). 

§2.  The  Hausdorff  metric  is  evidently  a distance  of  type  (in  the  space  of 

support  functions) . It  would  be  of  considerable  interest  to  find  an  appropriate 
formulation . 

§3.  Theorem  1 has  many  variants  along  traditional  lines.  For  instance,  uniform 
convergence  in  a subinterval  can  be  asserted  under  weaker  conditions.  Moreover,  point-wise 
convergence  rates  can  be  derived  for  each  of  the  component  s(p,F(t)). 

Theorem  2 suggests  that  linear  methods  may  not  be  natural  for  approximation  in 
T[K]  . C.  de  Boor  has  pointed  out  that  this  might  be  expected  since  the  formation  of 
all  point-wise  sums  in  constructing  is  not  mirrored  in  the  structure  of  a 

single,  fixed  non-convex  K (which  does  not  contain  all  connecting  line  segments) . 

§4.  The  extended  excursion  into  Korovkin  theory  for  families  of  functions  can  perhaps 
be  avoided  by  an  appeal  to  the  abstract  machinery  of  Banach  lattices  (see,  for  instance, 

[12,  esp.  V.2]).  We  have  not  seen  a clear  way  to  do  this,  and  in  any  case  the  quantita- 
tive formulation  given  may  be  of  particular  use.  Although  it  is  not  explicitly  given  in 
the  text,  the  following  bound  seems  best  possible 

H(TF,F)  < w_(u)  [H(F(0)  ) + 1]  + Y(F,T) 

— r 

where 

U2  = sup  ||  [T [ (t-x)  2B)  1 (x)  || 
x 

and 

y (F,T)  = sup(H(TG,G)  | G(x)  5 F(t)  , 0 < X < 1,  t 

Condition  (ii)  of  theorem  3 has  the  equivalent  (but  apparently  weaker) 
in  the  plane  (d=2)  : 


fixed)  . 
formulation 
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